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Exciton-polariton  band  structure  in  quantum-dot  lattices 

E.  L.  Ivchenko f,  Y.  Fu$  and  M.  Willanderf 
f Ioffe  Physico-Technical  Institute,  194021  St. Petersburg,  Russia 
J Physical  Electronics  and  Photonics,  Department  of  Microelectronics  and  Nano- 
science, University  of  Gothenberg  and  Chalmers  University  of  Technology, 
Fysikgrand  3,  S-412  96  Goteborg,  Sweden 


Abstract.  We  develop  a theory  of  exciton  polaritons  in  three-dimensional  quantum-dot  lattices 
with  the  period  comparable  to  the  light  wavelength.  A system  of  the  Maxwell  equations  and 
nonlocal  material  relation  are  used  to  derive  the  dispersion  equation  in  a rather  general  and 
well-converging  form.  A possibility  of  analytical  description  of  the  dispersion  is  questioned  and 
discussed.  The  photon  band  structure  is  calculated  for  a face-centered -cubic  lattice  with  spherical 
dots  of  the  radius  exceeding  the  bulk-exciton  Bohr  radius.  The  dispersion  along  the  E — X 
and  F — L lines  is  characterized  by  a strong  anticrossing  between  bare  transverse  photon  and 
exciton  branches  and  by  remarkable  overlapping  band  gaps.  Approaching  the  U and  W points 
the  exciton-polariton  branches  converge  and  the  gap  becomes  negligible. 


Recently  van  Coevorden  et  al.  [1]  have  calculated  the  optical  band  structure  of  a 
three-dimensional  (3D)  lattice  of  resonant  two-level  atoms.  They  have  solved  numeri- 
cally the  dispersion  equation  for  light  waves  in  a face-centered-cubic  atomic  lattice  and 
demonstrated  that,  in  the  certain  range  of  parameters,  there  exists  an  overlap  of  pho- 
tonic band  gaps  in  all  directions  in  the  frequency  region  near  the  two-level  resonance. 
Here  we  consider  the  photonic  (or  more  precisely,  exciton-polaritonic)  band  structure 
of  lattices  formed  by  a 3D  periodic  array  of  quantum  dots  (QDs). 

We  start  from  the  Maxwell  equations 

/U)\2 

AE  — grad  div  E = — ( — j D , 

div  D = 0 (1) 

for  the  electric  field  E and  the  displacement  vector  D.  The  nonlocal  material  equation 
relating  D and  E is  taken  in  the  form  (see  [2]) 

D(r)  = efeE(r)  + 47rPexc(r)  , (2) 

47rPm.(r)  = T(a,)^<f>a(r)  f <E>a(r')E(r>r  . (3) 

a 4 

Here  a are  the  lattice  translation  vectors  enumerating  quantum  dots,  <f>a(r)  = <f>0(r  - a) 
is  the  envelope  function  T't,j:c(r(,,  rh ; a)  of  an  exciton  excited  in  the  ath  QD  at  coinciding 
electron  and  hole  coordinates:  <f>a(r)  = T'an.(r,  r;a).  The  other  notations  are 

T(cu)=27T — 2 — > (4) 

— t oL 
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u>lt  and  aB  are  the  exciton  longitudinal-transverse  splitting  and  Bohr  radius  in  the 
corresponding  bulk  semiconductor,  u>Q  is  the  QD-exciton  resonance  frequency,  e*  is 
the  background  dielectric  constant  which  is  assumed  to  coincide  with  the  dielectric 
constant  of  the  barrier  material.  In  the  following  we  neglect  the  overlap  of  exciton 
envelope  functions  and  4/a,  with  a 7^  a'  so  that  excitons  excited  in  different  dots 
are  assumed  to  be  coupled  only  via  electromagnetic  field. 

It  follows  from  Eq.  (2)  that  div  E = —(47 r/e*)  div  Pexc  which  allows  to  rewrite  the 
first  Eq.  (1)  as 

AE  + k2 E = —47 (l  + k~2grad  div  ) Pexc  , (5) 

where  ko  — oj / c,  k = = c ont/c  and  = y/ei,. 

We  seek  for  Bloch-like  solutions  of  Eq.  (5)  satisfying  the  translational  symmetry 
EK(r  + a)  = exp(z'Ka)  EK(r)  where  the  wave  vector  K is  defined  within  the  first 
Brillouin  zone.  The  exciton-polariton  dispersion  w(K)  can  be  shown  to  satisfy  the 
equation 

Det||<5Q)3 -Ra0(uj,  K)||  = 0 , (6) 

where  a,  [3  = x,y,z,  5a,p  is  the  Kronecker  symbol  and,  for  QD  lattices, 


£qT(u;)  |fK+bl2^q/3(K  + b) 

vo  Y lK  + b<2  “ kl 


(7) 


Iq=  I Mr)eiQrdr  , SQ0  = Sa0  - , (8) 

b are  the  reciprocal  lattice  vectors  and  vo  is  the  volume  of  the  lattice  primitive  cell. 

Similarly  to  [1]  we  consider  a face-centered-cubic  lattice  with  the  lattice  constant 
a and  the  unit-cell  volume  vo  = a3/ 4.  It  is  convenient  to  introduce  a dimensionless 
parameter  P = (TrV^c/atuort/,)3  and  the  dimensionless  frequency  O = u>/u>q.  The 
calculation  is  performed  for  spherical  QDs  with  the  radius  R exceeding  the  Bohr  radius 
aB  in  which  case  we  have 


/ 2R\3/2  sin  QR 
Q_7rV ~b)  QRW  (: QR )2]  ' 

Then  Eq.  (7)  can  be  transformed  into 

Q2 

Ra0(n,  K)  = N^—j*a0(n,  K)  , 


(9) 


(10) 


cra0(Cl,  K) 


v/([K  + b|^)5ag(K  + b) 
Y 02-f72(K  + b) 


(11) 


N _ 64  Lou  fRV 

ir  too  \a ) 


. /w 


7r2  sinx  1 

x{i r2  — x 2)  ’ 


(12) 


fi(Q)  = cQ/ujQttb.  Eq.  (6)  is  equivalent  to  the  three  separate  equations  Rj(kl,  K)  = 1 
where  Rj  (/  = 1 , 2,  3)  are  eigenvalues  of  the  matrix  Ra0.  For  high-symmetry  points  of 
the  Brillouin  zone,  the  symmetry  imposes  certain  relations  between  the  Ra0  components 
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Table  1.  Dispersion  equations  written  in  terms  of  Rap  for  different  K points  in  the  Brillouin  zone. 


K (2-ir/a) 

Nonzero  components  of  Rni 

Dispersion  equations 

r (o,  o,  o) 

Rxx  — Ryy  — Rzz 

Rx. r=l 

A (0,0,1) 

Rxx — Ryy, 

/Tt=l,  Rzz=  1 

L (1/2,  1/2,  1/2) 

Raa=Rxx,  Ra(3=RXy(a'  7^  (3) 

Rxx  — Rxy=  1 , Rxx+2Rxy=  1 

IF  (1/2,  0,1) 

Rxxi  Ryy — Rzz 

Rxx=\,  ^yy=l 

K (3/4,  0,  3/4) 

Rxx — Rzz-,  Ryy  ? Rxz — Rzx 

RXX±RXZ=\,  Ryy=  1 

U (1/4,  1/4,1) 

Rxx — Ryy  , Rzz , Rxy — Ryx 

Rxx±Rxy=l,  R.:=  1 

Fig  1.  Exciton-polariton  dispersion  near  the  exciton  resonance  frequency  wo  in  a face-centered- 
cubic  lattice  of  spherical  QDs  characterized  by  the  following  set  of  parameters:  P = 1.1,  R/a  = 
1/4  and  wir/wo  = 5x  10~4.  The  dashed  lines  show  the  photon  dispersion  in  the  empty  lattice, 
i.e.  for  lolt  = 0,  the  dotted  horizontal  line  indicates  the  value  w = wo. 


and  the  eigenvalues  R,  can  be  readily  expressed  via  these  components  as  illustrated  in 
Table  1 for  the  points  F,X,L,  W,K  and  U. 

Further  simplifications  follow  taking  into  account  a small  value  of  the  parameter  N 
in  Eq.  (10)  since,  in  semiconductors,  the  ratio  lolt/lo0  typically  lies  between  10  4 and 

10  3.  Then  in  the  frequency  region  given  by  the  condition  |fi 1 P1/3  — 1 one  can 

readily  use  the  approximate  equation  H-1k  (N /2)oy(l,  K)  where  Uj  are  eigenvalues 
of  the  cra0  matrix. 

Fig.  1 shows  the  photonic  band  structure  calculated  for  the  dots  of  radius  R = a/ 4 
and  forP  = 1.1,  ujlt/ujq  = 5 x 10“4.  The  dispersion  on  the  A line  is  characterized  by 
a giant  anticrossing  between  the  branches  of  bare  transverse  photon  and  exciton  modes. 
At  the  X point,  the  gap  is  determined  by  the  separation  between  the  longitudinal  and 
lower  transverse  branches,  it  is  still  remarkable  and  exceeds  0.5uy.7\  However  near  the 
points  U and  W the  excitron-polariton  branches  converge  and  the  gap  almost  disappears. 
Note  that  the  anticrossing  can  be  described  with  a high  accuracy  by  retaining  in  the 
sum  over  b in  Eq.  (11)  the  two  terms  due  to  b = 0,  -(4-7r/a)(0,  0,  1)  for  the  A points 
and  b = 0,  -(27r/a)(l,  1,  1)  for  the  A points. 

The  3D  QD  arrays  with  periods  comparable  with  the  light  wavelength  (P  « 1) 
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and  with  sizes  exceeding  the  bulk-exciton  Bohr  radius  could  he  grown  artificially  or 
by  embedding  semiconductor  microcrystals  into  the  pores  of  porous  materials  like  the 
synthetic  opal  [3].  It  should  be  mentioned  that  the  developed  theory  takes  into  ac- 
count a contribution  of  only  one  exciton  resonance  which  is  valid  if  the  separation 
between  the  exciton  size-quantization  levels  is  much  larger  than  the  bulk  value  of 
the  exciton  longitudinal-transverse  splitting,  wLT.  In  the  opposite  limit  of  extremely 
large  bulk-exciton  translational  effective  mass  one  can  use  the  local  material  relation 
D(r)  = e(r,  cu)E(r)  as  it  was  done  by  Sigalas  et  al.  [4]  for  phonon-polaritons  in  a 
two-dimensional  lattice  consisting  of  semiconductor  cylinders. 

E.  L.  I.  acknowledges  partial  support  from  Russian  Foundation  of  Basic  Research. 
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